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A Model-Centric Method And Apparatus For Dynamic 
Simulation, Estimation And Optimization 

1 . Field of the Invention 

This invention relates to dynamic process modeling 
and more particularly to dynamic simulation, estimation 
and optimization of industrial processes with a flexible 
model development framework. 

2 . Description of the Prior Art 

Facing economic, environmental and safety 
10 constraints, industries are forced to improve performance 
of their plants and to reduce the cost and time for plant 
development. As reported by Marquardt, W. (1991), 
"Dynamic process simulation - recent trends and future 
challenges" in Arkun, Y. and W. H. Ray, (Eds.) Chemical 
Process Control CPC-IV CAChE, Austin, AIChE, New York, 
pl31-180; and Pantelides, C. C. and P. I. Barton (1993),, 
"Equation oriented dynamic simulation: Current status and 
future perspectives", Computer and Chemical Engineering, 
17S, p263-285, system simulation, estimation and 
20 optimization have become vital tools in order to achieve 
these goals and to be competitive in a rapidly growing 
market with increasingly short innovation cycles. 

For example, consider the pulp and paper industry 
where continuous pulping processes have, as is described 
by Gullichen, J., and C.-J. Fogelboln (2000), Chemical 
Pulping, Gummerus Printing, Jyraeskylae, Finland, grown 
to dominate the production of kraft pulp. The digester 
is the heart in a pulp and paper mill, and its operation 
directly impacts the economic benefit of the plant and 
30 the quality of the final paper product. Therefore, the 
digester has attracted much attention in model 
development and advanced control etc. due to its critical 
importance and high capital investment. See Wisnewski, 
P- A., F. J. Doyle III, and F. Kayihan (1997), 
"Fundamental continuous pulp digester model for 
simulation and control", AIChE J., 43(12), p3175-3192; 
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and Dahlquist, E. et al.(2001), "Economic benefits of 
advanced digester control", Pulp Digester Modeling and 
Control Workshop, June 28, 2001, Annapolis, MD, USA. 

While simulation techniques have been commonly used 
in the process industries such as the petrochemical 
industry for several decades, there still are many 
challenging issues to be solved. In particular, the pulp 
and paper industry has had much less applications of this 
technology to date. A few steady-state process 

10 simulation packages are available for the pulp and paper 
industry, but these are sequential modular simulators and 
are customized for specific needs in the industry. 

Although some extensions to dynamic simulation of 
the pulping process using one of these packages have been 
made [see Shewchuk, C. F. (1987), "Massbal mkll: New 
process simulation system", Pulp & Paper Canada, 88(5), 
p76-82], the pulp and paper industry still lacks a 
dynamic modeling system with rigorous process modules, 
especially for a digester. On the other hand, interest is 
20 increasing in efficient dynamic simulation systems 
integrated with advanced application tools like 
simulation, estimation and optimization for that 
industry. 

In any chemical process industry, there is a constant 
drive to maximize throughput, minimize product quality 
fluctuations, etc. amidst constant demand for varied 
product grades. In such a situation increased 
understanding of the underlying process and enhanced 
knowledge of the operations becomes extremely important. 
30 Unfortunately, due to the nature of the operations and 
equipment, it is not always possible to access all 
sections of the process vessels. 

For example, in the digester where the cellulose 
fibers are separated from lignin, the wood chips are 
steamed and fed with cooking chemicals in the form of 
alkali with very long residence times (about 6 to 8 
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hours) . It is important to know the amount and extent of 
reaction inside the digester at various points, so that 
the side draws can be manipulated to attain the 
operational targets. However it is very difficult or even 
impossible to obtain accurate measurements inside the 
digester. The current practice is to take samples at the 
side draws and analyze them with online analyzers or lab 
experiments. Lab results are time consuming and not 
available frequently. Online analyzers based on side draws 

10 provide continuous measurements but provide only local 
measurements and do not provide information inside the 
digesters. Moreover they are susceptible to error or even 
failure and may require high maintenance. 

In general, many engineering applications such as 
process design, simulation, estimation, and optimization 
in the chemical process industries have been implemented 
separately with different software. These applications 
rely on the presence of a good mathematical model, but 
the models used may be developed from different sources 

20 in spite of the same process. These approaches result in 
increased project engineering costs due to redundant 
effort as there is little or no model reuse and design 
changes are not easily translated across all models. 

The success of the model-based applications depends 
largely on the approximation accuracy of the dynamic 
model. Usually a rigorous model is preferred, which must 
represent complex physical phenomena within this process 
including both chemical reactions and transport phenomena 
such as bulk flow, convection and diffusion. However 

30 these models are often very large and complex, involving 
thousands of nonlinear algebraic and differential 
equations. For example, in the pulping process described 
above, the digester is a tubular multiphase vertical 
reactor operating in plug flow with varying directional 
flows of liquor throughout the vessel, and presents quite 
some challenges for process model development. Various 
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configurations of digesters also exist within the 
industry. 

It would be desirable for a flowsheet simulator to 
have the ability to handle these different types and 
geometries of operation. This also means that a flexible 
and efficient model development framework for simulation 
and other model-based solutions is required. The present 
invention meets these requirements as it presents a 
framework that allows a process flowsheet to be built for 

10 dynamic simulation, estimation and optimization based on 
a model-centric methodology. 

Summary of the Invention 
A method associated with a process. The method 
comprises building a general purpose flowsheet for the 
process. The general purpose flowsheet representing a 
model of the process. The method also comprises 
generating the model of the process from the general, 
purpose flowsheet. The method further comprises 

executing upon the model generated from the general 

20 purpose flowsheet applications selected from simulation, 
estimation and optimization. 

A computer readable medium having instructions 
associated with a process. The instructions comprise 
building a general purpose flowsheet for the process. 
The general purpose flowsheet representing a model of the 
process. The method also comprises generating the model 
of the process from the general purpose flowsheet. The 
method further comprises executing applications selected 
from simulation, estimation or optimization upon the 

30 model generated from the general purpose flowsheet. 

A method for enabling a user to build a general 
purpose process flowsheet whose underlying model may be 
used to execute applications selected from simulation, 
estimation and optimization. The method comprises 

enabling the user to build the general purpose process 
flowsheet by providing a graphical flowsheet builder; a 
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library of common equipment models; selectable chemical 
species and their thermophysical property definitions; 
and stream connectors to allow the user to connect models 
selected from the library. The method also comprises 
providing an interpreter to generate the underlying model 
from the general purpose flowsheet. The method further 
comprises providing a solution engine that enables the 
user to use the generated underlying model in 
applications selected from simulation, estimation and 

10 optimization. 

A method for enabling a user to build a general 
purpose process flowsheet whose underlying model may be 
used to execute applications selected from simulation, 
estimation and optimization. The method comprises 

enabling the user to build the general purpose process 
flowsheet by providing a graphical flowsheet builder; a 
library of common equipment models; selectable chemical 
species and their thermophysical property definitions; 
and stream connectors to enable the user to connect 

20 models selected from the library. The method also 
comprises providing an interpreter to generate the 
underlying model from the general purpose flowsheet. The 
method further comprises providing a solution engine that 
allows applications selected from simulation, estimation 
or optimization to be executed upon the generated 
underlying model. 

An apparatus for building a general purpose 
flowsheet having an underlying model for a process. The 
apparatus comprises a library of selectable common 

30 equipment models. The apparatus also comprises a 

repository of selectable chemical species and their 
thermophysical property definitions. The apparatus 

further comprises a flowsheet editor for selecting one or 
more of the selectable common equipment models and one or 
more of the selectable chemical species and 
thermophysical properties definitions associated with the 
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one or more selected models and linking the one or more 
selected models by stream connectors to build the general 
purpose flowsheet. The apparatus also further comprises 
an interpreter to generate the underlying model from the 
general purpose flowsheet; and a solution engine for 
executing applications selected from simulation, 
estimation or optimization upon the generated underlying 
model . 

Description of the Drawing 
10 Fig. 1 shows a simplified schematic of an industrial 

pulping process. 

Fig. 2 shows the configuration architecture for the 
modeling framework or model development framework 
described herein. 

Fig. 3 shows the execution architecture for the 
modeling framework or model development framework 
described herein. 

Fig. A shows a screen shot for and Fig. 4a shows in 
detail the equivalent representation of the pulp digester 
20 as a flowsheet in the model development system and 
framework described in this embodiment. 

Fig. 5 shows in an investigation of the dynamic 
behavior of the model developed in accordance with the 
present invention the movement of the Kappa number profile 
after a predetermined step increase of the chip screw 
while keeping all other process variables unchanged. 

Fig. 6 shows a comparison of the bottom Kappa number 
between the prediction of the model developed in 
accordance with the present invention and measurement over 
30 a day. 

Fig. 7 shows based on the model developed in 
accordance with the present invention the prediction of 
the moisture in the chip feed entering the digester. 
Description of the Preferred Embodiment (s) 

The present invention is described below using the 
pulp digester, a simplified schematic for which is shown 
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in Fig. 1, as an application case. It should, however, be 
appreciated that the model-centric methodology of the 
present invention may be used to dynamically simulate, 
estimate and optimize other processes. 

As is shown in Fig. 1, and as is well known to those 
in the art, in the pulping process 10 wood chips are mixed 
with steam in the chip bin 12 and are cooked in the 
digester 14 to separate the lignin from the cellulose 
fibers to form pulp. Digester 14 can in accordance with 

10 the present invention and as is described in more detail 
below be represented using four process modules, namely, 
top separator 14a, digester sections 14b, screens 14c and 
bottom scraper 14d. 

The present invention is described in connection with 
a modeling framework (or. a model development framework), 
the configuration architecture for which is shown in Fig. 
2 and the execution architecture for which is shown in 
Fig. 3, that supports graphical user interfaces enabling 
component configuration and editing, graphical process 

20 f lowsheeting, model specifications, steady state/dynamic 
simulation configuration, and integration with estimation 
and optimization applications. 

The framework 20 of Figs. 2 and 3 allows a user to 
build a process flowsheet, such as the flowsheet shown in 
Figs. 4 and 4a which shows the digester section model 50 
of pulping process 10, for simulation of industrial 
processes, and use this same flowsheet in a variety of 
applications such as inferential property calculations, 
parameter estimation, data reconciliation, optimization, 

30 etc. It provides an easy to use graphical flowsheet 
builder (editor) , a library of common equipment models, 
and flexible interfaces for configuring activities or 
applications . 

As is shown in the configuration architecture of Fig. 
2, a process flowsheet such as the flowsheet shown in 
Figs. 4 and 4a, is built using the 'flowsheet editor' 22 
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by dragging and dropping equipment shapes in the form of 
icons in flowsheet editor stencils representing 
mathematical models from the model library 24 and 
connecting these shapes by stream connectors. Streams are 
composed of the different chemical components configured 
using the 'component editor' 26. Editor 26 can link to a 
framework-provided properties databank 21a or has an 
interface to other commercially available property 
packages such as DIPPR (21b) to obtain chemical component 

10 property information. The component editor 26 is also used 
to specify the thermodynamic package 41 (see Fig. 3) used 
for calculating thermodynamic properties. Simulation, 
estimation, and optimization activities are configured 
using the corresponding 'activity editors' 28, 30 and 32, 
respectively. All configuration information using the 
various editors are persisted in the Flowsheet Databases 
23 of Figs. 2 and 3. 

The task scheduler 43 of Fig. 3 commands the 
estimation, optimization and simulation 'input generators' 

20 38, 40, 42, respectively of Fig. 3, to collect all 
necessary information from the persisted storage 
(flowsheet database and model library files) and submit 
them to the solution engine 44 for solving the resulting 
system of flowsheet models and equations. When operating 
data is required for example in estimation activities 
either to tune the models or to update the model 
parameters periodically, an interface to the historian 
available with the framework can link to the plant 
historian 45 to fetch required data. Such data 

30 connectivity for model variables from various data sources 
is provided by the 'data source editor' 36 of Fig. 2 
allowing completed models to be parameterized and 
validated against available data. The 'variable editor' 
25 of Fig. 2 provides an interface into all of the 
equipments, variables, and parameters to show the values, 
engineering units, and upper and lower bounds on the 
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variables . 

The framework enables the use of the validated models 
in a number of off-line and on-line applications including 
"what-if" analyses, engineering design, real-time 
simulation, inferential property prediction, data 
reconciliation, yield accounting, soft-sensor analysis, 
control applications and process optimization. The main 
benefit of the framework is its model-centric capability 
that allows a single process model to be deployed in a 

10 multiple number of applications. 

The equipment of the flowsheet are represented using 
objects in the underlying database. These objects also 
provide their data via an OPC server to external systems 
and are also integrated into a suitable control platform 
such as the Industrial 11 architecture of the Aspect 
Integration Platform available as of the filing date of 
this U.S. patent application from the assignee of the 
present invention. Once the flowsheet configuration is 
complete in the framework, the control platform objects 

20 can be created simply using an uploader tool. Conversely, 
the flowsheet may also be created in the framework from 
the collection of selected objects in a control system 
using a simple tool. 

The model development framework is flexible to 
support different solution engines such as sequential 
modular solvers and equation-oriented solvers. By 
default, it includes a solution engine 44 that is an 
equation-oriented simultaneous solver for differential- 
algebraic equation (DAE) systems. In the embodiment for 

30 the present invention described herein this solution 
engine is the gPROMS software, which is available as of 
the filing date of this U.S. patent application from 
Process Systems Enterprise Limited of London, England. 

The results of executing various activities in the 
framework can be obtained in a variety of tools including 
some provided by Process Systems Enterprise such as gRMS, 
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45a (a graphical results management system) , gExcelOutput 
45b (results sent to spreadsheet), and gPlot 45c (results 
stored in files) . The framework provides additional 
methods to view the results of activities. Preconf igured 
Preformatted reports 4 6 (for example, equipment reports 
and measurement reports) can automatically be generated 
in Microsoft® Excel at the end of an execution. The 
results are also available via an inbuilt OPC server 44a 
in solution engine 44 to external systems. 

10 The combination of the model development framework 

and the solution engine gives rise to a state of the art 
general purpose process flowsheet modeling and simulation 
environment. This framework is used to develop a 
rigorous dynamic model of the digester 14. Then based on 
the same model, the model-based applications are 
implemented. For example, periodic simulation' and 
estimation using the model can infer the difficult to 
measure variables and unmeasured disturbances that enter 
the pulping process 10. The framework allows for 

20 extensive model reuse and simple change propagation to 
all applications and thus reduces the total project 
engineering effort. 

In order to develop a dynamic process model 
efficiently and flexibly, a model library 24 (see Fig. 2) 
is firstly developed within the framework. For the kraft 
pulping digester process 10, a unified modeling approach 
for the pulp digester 14 is presented through a 
modularization concept to meet different configuration 
needs and enhance reusability of the model. 

30 As described in connection with Fig. 1, four process 

unit modules, Digester Sections 14b, Screens 14c, Top 
Separator 14a and Bottom Scraper 14d are considered for 
the digester 14. The digester model can easily be 
configured by connecting these basic unit modules. The 
Digester Section module 14b is formatted based upon a 
nonlinear DAE model by lumping the distributed model. 
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This allows for various levels of complexity balanced by 
efficiency (accuracy and size) of solution through 
selecting the number of the Digester Sections. 

All process modules in the model library 24 are 
developed based upon the first-principles approach. 
Rigorous reaction kinetics, dynamic heat and mass 
balances, momentum transport equations and material and 
thermal dif f usivities are all considered in the model 
equations. Thermodynamic packages are linked within the 

10 framework to provide physical property information to the 
model components. 

Each block or process module on the flowsheet of 
Figs. 4 and 4a represents an equipment or a facility such 
as a valve, header, or a pump. The lines connecting these 
blocks are the streams carrying component information 
such as composition, pressure, temperature, and flowrate. 
The unconnected blocks represent facilities to include 
custom defined calculations that can be used to perform 
user desired calculations outside the scope of the 

20 standard library models. 

With the aid of the models and the framework 
described above, various steady state and dynamic 
simulations can easily be configured. Specifications can 
be configured based upon the physical measurements in a 
mill, design data corresponding to the equipment 
dimensions, and some empirical data. Through the 

modeling framework, the dynamic solution of the rigorous 
process model can be solved quickly and robustly. The 
process simulation allows various investigations 

30 supporting both process studies and engineering analysis. 

To investigate the dynamic behavior of the model, 
some tests were made through making step changes of 
several process variables. Figure 5 shows that the Kappa 
number profile moves after making a 25% step increase in 
the rotations per minute (rpm) of the Chip screw, while 
keeping other variables unchanged. The Kappa numbers at 
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all locations along the vessel increase with different 
gains and rise times. 

The developed digester process model needs to be 
tuned and validated against operation data. This is 
important for two reasons. First, there are a number of 
parameters, such as kinetics parameters and heat transfer 
coefficients, which are unique to a particular process. 
Second, any model should be well proven before it gains 
the political acceptance required to support major 
engineering decisions in an operating plant. In 
accordance with the present invention the model 
parameters are updated by estimation based upon the 
collected operation data. 

A comparison of the bottom Kappa number between 
model predication and measurement over a day are shown in 
Figure 6. The figure shows that the model prediction 
results can reasonably match with actual plant 
measurements . 

Before performing model-based applications, dynamic 
data reconciliation is used to reconcile process and 
laboratory measurements from the plant to satisfy 
material, energy and momentum balances. The parameter 
estimation described above can also be performed 
simultaneously with data reconciliation. The framework 
performs a rigorous calculation to solve the estimation 
and reconciliation problem and as is described above uses 
a solution engine such as the gPROMS software. 

The parameter estimation and data reconciliation is 
setup as an optimization problem in which an attempt is 
made to find the values of all unknown estimated 
parameters (including measurement biases) such that the 
mathematical model will predict the values obtained from 
the plant operation. A maximum likelihood objective 
function is used to perform parameter estimation and data 
reconciliation. The configuration of the reconciliation 
problem can be carried out through the graphical 
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estimation interface 30 of Figure 2. Based on the well- 
tuned model, it is possible to infer the difficult to 
measure variables and unmeasured disturbances that enter 
the process. Figure 7 shows the prediction of the 
moisture in the chip feed entering the digester. 

While the framework of the present invention is 
described above with respect to a continuous process it 
should be appreciated that the same framework may also be 
used for batch or other types of processes. 

It is to be understood that the description of the 
preferred embodiment (s) is (are) intended to be only 
illustrative, rather than exhaustive, of the present 
invention. Those of ordinary skill will be able to make 
certain additions, deletions, and/or modifications to the 
embodiment (s) of the disclosed subject matter without 
departing from the spirit of the invention or its scope, 
as defined by the appended claims. 



